Quantum coherence control of solid-state charge qubit by means of a suboptimal 

feedback algorithm 
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The quantum coherence control of a solid-state charge qubit is studied by using a suboptimal con- 
tinuous feedback algorithm within the Bayesian feedback scheme. For the coherent Rabi oscillation, 
the present algorithm suggests a simple bang-bang control protocol, in which the control parameter 
is modulated between two values. For the coherence protection of idle state, the present approach 
is applicable to arbitrary states, including those lying on the equator of the Bloch sphere which are 
out of control in the previous Markovian feedback scheme. 
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Introduction. — Being stimulated by the interest of 
solid-state quantum computation, in recent years the 
measurement problem of solid-state qubit is becoming 
an intensive research subject. In particular, the setup of 
a charge qubit measured by a mesoscopic transport de- 
vice, for instance, the quantum point contact (QPC), has 
received considerable attention 0, S 0, 13 • Very recently, 
the mesoscopic QPC was experimentally demonstrated 
as a charge-sensitive electro- meter 

Strikingly, rather than the projective strong measure- 
ment, the continuous weak measurement can be em- 
ployed to realize some marvellous tasks, such as quantum 
state initialization, stabilization, and coherence protec- 
tion, etc IE ■ The central idea is the measurement- 
result-based feedback control, which falls typically into 
two main categories, say, the Markovian scheme |3] and 
the Bayesian scheme |2, E|. In quantum optics, the 
study of quantum feedback control has been going on 
for more than a decade. However, in solid states it is 

a relatively new subject mill mill- 

In this work, 

based on a suboptimal feedback algorithm 01 > study 
the problem of quantum coherence protection of a solid- 
state charge qubit which is subject to the continuous 
weak measurement by a mesoscopic QPC. 

Formalism. — To be specific, the charge qubit is as- 
sumed to be two coupled quantum dots (CQD's), occu- 
pied by a single electron [l|, |l^ [13 , whose Hamiltonian 
reads = e(cjc2 — cjci) -I- Vl{c\c2 + c\ci). Here cj (q) 
is the electron creation (annihilation) operator in the ith 
dot, £ and are, respectively, the energy offset of the 
two local dot-states and their coupling strength. At zero- 
temperature and in the limit of large bias voltage across 
the QPC, the measurement back-action on the qubit is 
described by the following master equation {fi = 1) il3] 



p{t) = ^i[Hqb, p{t)] + V[T + Xni]p{t) 



(1) 



where ni = c|ci is the occupation- number operator of 
the first dot, the parameters T and X are given by |Tp = 
2T^\T\^gigreV = D2, and \T + X]-" = 2n\T + x?gi9reV = 
Di. Here, the tunneling amplitudes through QPC, T 



and Xj 3're assumed to be independent of the reservoir 
states, and real for simplicity; gi and gr are the density of 
states for the left and right reservoirs; V — [pL — II'r)/^ 
is the bias voltage across the QPC. The superoperator 
in Eq. (Q) is defined as V[r]p = J[r\p — A[r]p, where 
J[r\p = rpr\ and ^M/O = (r^rp -|- pr'^r)/2. Intu- 
itively, corresponding to the qubit electron in either the 
first or the second dot, the QPC measurement current 
of Ji = eDi or I2 = eD2 is expected. However, for 
continuous weak measurement, owing to the qubit elec- 
tron jumps between the two dots, the continuous output 
current I[t) is a complicated stochastic random quan- 
tity. But, quite remarkably, by using I(t) one can de- 
sign a proper feedback Hamiltonian Hfbit) to maintain 
or improve the quantum coherence of the qubit. Condi- 
tioned on Ijt) and involving Hfh(t), Eq. Q is unravelled 

as miiiiia 



dpc{t) 



-i[Hqb:Pc]dt + ^ "fjV[(7.j\pcdt 



+V[§a,]p,dt + H[^a,]pdW{t) 
'~i[Hfb{t),Pc]dt. 



(2) 



Compared to Eq. , in addition to the feedback term 
(the fifth one), we have included here the possible noise 
term (the second one) due to external environment, 
which would spoil the quantum coherence and needs 
to be eliminated by the feedback control strategy. In 
Eq. ||2J), the superoperator "7i" is defined by 7i[r]p — 
rp + pr^ — p Tr[rp + pr^. dW{t) is the Weiner increment, 
which is a Gaussian-white-noise stochastic variable, hav- 
ing the property of E[dW{t)] = and E[dWit)dW{T)] = 
S{t — T)dt, with E[- ■ ■ ] the ensemble average over a large 
number of stochastic realizations. The Weiner increment 
dW{t) is associated with the output current I{t) in terms 
of the relation: 



I{t) - I = Alipl, - p^2)/2 + ./S^m, (3) 



where (,{t)dt - 
and 5*0 = 2el. 



dW{t), 



/ = (/i +/2)/2, A/ = /i - I2, 
{i\pc{t)\j) is the density matrix 
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FIG. 1: The ensemble evolution of pii{t) = E[p1i{t)] for dif- 
ferent feedback strengths /i=0 (gray solid curve); 0.001 (black 
solid curve); 0.1 (dashed curve). 



element over the qubit basis states. Based on the out- 
put current I{t), there are two typical schemes to design 
the feedback Hamiltonian, i.e., the so-called Markovian 
"/(t)-based" feedback 6], and the Bayesian "pc{t) -based" 
feedback ^ J^] . In the following, the second one will be 
implemented to protect the qubit coherent evolution and 
its idle state, on the basis of a suboptimal feedback algo- 
rithm 15]. 

Feedback Control of Rabi Oscillation. — For simplicity, 
consider the symmetric qubit (e = 0). The "desired" 
pure-state evolution to be protected is 



\Mt)) - cos(m)|l) 



i{nt)\2) 



(4) 



where |1) and |2) are the dot states. The basic idea of 
Bayesian feedback control is to carry out first an esti- 
mate of the qubit state pc(t) based on the noisy current 
I{t), then compare it with the desired state ipdit) and 
design the feedback Hamiltonian using the calculated dif- 
ference. In Ref. ^3 the estimated quantity is the relative 
phase between the superposed states "|1) and |2)", i.e., 
(j){t) = arctan(2Imp^2(0/[Pii(0 ~P22(*)])> which is com- 
pared with (j}Q = 2nt{ mod 27r) defined from the desired 
state. Then, the phase "error" A(f){t) — (f){t) — (f)Q is used 
to design the feedback Hamiltonian, i.e., Hfb{t) oc A(/)(t). 
In certain sense, this scheme may be inconvenient to be 
implemented in practice, since the feedback Hamiltonian 
is a complicated stochastic function of time. In the fol- 
lowing, we propose an alternative scheme based on a sub- 
optimal algorithm [T^ . 

For real-time feedback control, each feedback step acts 
only for an infinitesimal time, "A<" . The suboptimal 
algorithm is based on the principle that the state evo- 
lution in each infinitesimal time step will maximize the 
fidelity of the estimated state with the desired (target) 
state. The state evolution in the presence of feedback 
is governed by Eq. ||2J). As far as the term related to 
the feedback Hamiltonian is concerned, the final state 
Pc{t + At) is given by 



Pc{t + M) ^ pc{t) -t[Hfb,pc{t)]At 

-i[i7/,,[i//b,p,(0]](Ai)2 



FIG. 2: Synchronization degree D as a function of feedback 
strength p, for several magnitudes of dephasing due to envi- 
ronment: 7z = (solid curve); O.OSF^ (dashed curve); and 
0.25F£j (dash-dotted curve). 



The fidelity of this state with the target state is defined 
by 

{i;d\pcit + At)\'tPd) = {i^d\pc{tMd) 

- t{iPd\[Hfb,Pc{t)]\iPd)At 

- ^{U[Hfi„[Hfb,Pc{t)md){At)^ 
+ ■•• . (6) 

To optimize the fidelity, one should maximize the coef- 
ficient of At, which is the dominant term. Similar to 
other control theories, the maximization must be sub- 
ject to certain constraints, such as the restriction on the 
maximum eigenvalue of -ff/b, the sum of the norms of 
the eigenvalues, or the sum of the squares of the eigen- 
values, etc. Physically, these constraints stem from the 
limitation of the feedback strength or finite Hamiltonian 
resources. Here we adopt the last type of constraint, 
namely, 



(7) 



Under this constraint, the feedback Hamiltonian can be 
constructed as UM 



Hfb^ic[\Mt)){Mt)lPc{t)], 



(8) 



with a ^ ('ijjdit)\pl{t)\xjjd{t)), and 



where c =- y ^(^^ . 

6 = {ipd{t)\pc{t)\i)dit)). Substituting Eq. (g)) into Eq. 
yields 



with 



^/f , A<^ < 



(9) 



(10) 



(5) 



In this case, the suboptimal feedback control resembles a 
simple bang-bang control, where the feedback parameter 
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FIG. 3: (a) 3D-plot of the synchronization degree D of the 
feedback control for the measurement ; (b) 2D-plot of D of 
the feedback control for several magnitudes of dephasing due 
to environment: 7z = (solid curve); COSFd (dashed curve); 
and 0.25rti (dash-dotted curve). The feedback strength is 
= 0.1. 



between two values is determined by the phase error A0 
as defined above. However, for other control purpose such 
as the protection of an arbitrary idle state which is to be 
studied soon, the suboptimal feedback scheme does not 
fall into the category of the "bang-bang" control, where 
the feedback Hamiltonian is more complicated. 

Now we present a Monte Carlo simulation for the feed- 
back effect based on the above algorithm. In the simu- 
lation the coupling strength between the qubit and the 
detector is characterized by the measurement-induced de- 
phasing rate, which is assumed Yd = X"^ /2 — 0.2557. 
In terms of the ensemble-averaged quantity pii{t) — 
E[p'ii{t)], Figure 1 shows that the simple feedback de- 
signed above can really eliminate the back-action of the 
detector, where the ideal Rabi oscillation is eventually 
recovered from the damped one by increasing the feed- 
back strength. The result in Fig. 1 is carried out for the 
case of ideal detector and in the absence of external noisy 
environment, i.e., by setting 7j — 0. 

To quantitatively characterize the effect of feedback, 
the so-called synchronization degree is introduced as , 
D = 2{TrpcPd) — 1. For instance, the complete syn- 
chronization is characterized by I? = 1, indicating that 
the desired state is perfectly protected by the feedback. 
In Fig. 2 we plot the synchronization degree versus the 
feedback strength /x, in the presence of external "CTz"- 
dephasing with rates 7^ = 0, O.OSF^j, and 0.25Fd, respec- 
tively. In certain sense, the "tTz"-type error caused by 
the external environment can be regarded as the con- 
sequence of a non-ideal detector since the observable is 
also az- Figures 1 and 2 show that the fidelity of such 
feedback can be arbitrarily close to 100%, in the case of 
using ideal detector and in the absence of external noises. 
Nevertheless, it decreases in the case of a non-ideal de- 
tector and/or in the presence of significant interaction 
with external environment, as well as in the case of finite 
bandwidth of the line which carries signal from the de- 
tector. This limitation also existed in the scheme of Ref. 



EH 



of using the above suboptimal feedback to protect the idle 
state of a qubit. This is essential in information storage, 
while the Rabi oscillation discussed above corresponds 
to information processing. In general, an arbitrary pure 
state of qubit can be expressed as 



IV't) =cos-|l)+sin-e'*|2) . 



(11) 



The goal of feedback is to make it immune against the 
influence of environmental noise. 

Following the same procedure used above, say, max- 
imizing the "At" term in Eq. ((H)), and substituting 
Eq. (jni) into Eq. ©, we obtain 



Hfb = \[ciaz + C2<yx + C^CTy] , 



(12) 



with 



ci = 2a/3(Repi2 sin </> + Ini/0i2 cos (j)) , 

C2 = (/3^ - a'^)lmpl2 - a(3{pl^ - P22) sin (/) , 

C3 = al3 cos(t){pl^ - p^2) - ("^ - /3^)Rep^2 > 

in which a = cos |, and /3 = sin |. 

So far the feedback Hamiltonian is constructed on the 
basis of maximizing the coefficient of At in Eq. ^ . How- 
ever, this algorithm does not work if the target state is 
the eigenstate of pc, namely, Pc\iPt) — ^t\iPt)- For in- 
stance, such situation occurs when 9 = and tt. In this 
case, the dominant "At" -term in Eq. © vanishes, and 
one is forced to maximize the coefficient of (Ai)^. Un- 
der the constraint on the Hamiltonian as described by 
Eq 10, the explicit construction of the optimal Hft, reads 

m 



Hfb = ^/JI/2{\vi){iPt\ + |V't)(^^i|), 



(13) 



where \vi) is the eigenvector with the largest eigenvalue 
of Pc{t). Notice that for the idle target state |1) (or |2)), 
the conditional evolution described by Eq. Q will not 
generate quantum coherence between |1) and |2), because 
there is no coherent driving source. Thus the states |1) 
and 1 2) must be the eigenstates of Pc{t), with eigenval- 
ues Pii(t) and P22(*); respectively. From Eo. (fT^ . the 
concrete form of feedback Hamiltonian is detailed as 



H 



fb 







for 61 = ( IV't) = |1)), and 



H 



fb 







Pll < P22 
Pll > P22 



Pll < P22 
Pll > P22 



(14) 



(15) 



Protection of Idle State. — Now we turn to the subject 



for 6 = n (Ii/it) = |2)). Here we have used the property 
that if \vi) in Eq. (|13|) is the target state, we should set 
Hfb = for that time step, since in this case the fidelity 
will not increase under the feedback action. 

In the following numerical simulations, we first con- 
sider the problem of merely eliminating the back-action 
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FIG. 4: 3D-plot of the synchronization degree D of the feed- 
back control for the Ux-tjpe error. The parameters are feed- 
back strength = 0.1 and noise strength ')x = COSF^. 

of an ideal detector (i.e. 7j = 0), in the absence of any 
other external influence of environment. The result is 
shown in Fig. 3(a). It is observed that the synchroniza- 
tion degree is independent of the relative phase "0" be- 
tween the two superposed components [c.f. Eq. Hll|) ]. but 
depends on "6'" symmetrically around "6* = 7r/2" (the 
minimal point). Since observing "ctz" will induce de- 
phasing of the target state, this leads to the smaller syn- 
chronization degree for the more coherently superposed 
state. For = and tt, the corresponding target states 
are the eigenstates of the observable cr^; they are thus 
immune against the back-action of the detector, leading 
to D = 1 and Hfi, = according to the proposed feed- 
back scheme. On the other hand, for = 7r/2, the cor- 
responding target state is |^t) = '^^\^) + 6"^|2)), which 
is the most coherently superposed state. This state is 
destroyed most seriously by the detector's back-action, 
whose synchronization degree is thus minimal. 

Now we consider the influence of external environment. 
In reality, a possible error is the az-type, which causes 
similar pure dephasing as the detector's back-action does. 
As a consequence, the "(0, (/))"-dependence behavior of 
the synchronization degree shown in Fig. 3 (b) is similar 



to Fig. 3(a), but it decreases with the dephasing strength. 
Another possible error is the CTj;-type, which causes relax- 
ation between the states |1} and |2}. The feedback result 
is shown by Fig. 4, where the synchronization degree de- 
pends on both "6'" and "(/>" . In particular, in contrast to 
the "iTz"-type error, here for the states Ii/it) = |1) and |2) 
(corresponding to 6* = and tt) the feedback efficiency is 
worst, since they are the most seriously affected states by 
the "cra;"-type error. Note also that the QPC tunneling 
amplitude related to |1) is smaller than that related to 
1 2), which leads to fe^o > Dg^y^. 

Conclusion. — To summarize, based on a suboptimal 
feedback algorithm we have presented a study on quan- 
tum coherence control of a solid-state charge qubit, for 
both its quantum evolution and idle state. For the co- 
herent Rabi oscillation, our study leads to a simple bang- 
bang control feedback scheme, in terms of modulating a 
single control parameter by only two values, which dif- 
fers from the existing Bayesian scheme and may be 
easier to implement in some cases. For the coherence pro- 
tection of idle state, the present approach is applicable 
to arbitrary states, including those lying on the equator 
of the Bloch sphere which are out oj control in the pre- 
vious Markovian feedback scheme 20]. In addition to 
eliminating the measurement back-action, the proposed 
feedback scheme can also protect the qubit from environ- 
mental influence to some extent. However, as many other 
approaches, it cannot protect the qubit from drastic in- 
fluence of external environment. This drawback may be 
partially overcome by developing other optimal feedback 
algorithms, which is a promising but challenging task in 
this field. 
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